Growing Ising-like chain as a model of online emotional interactions 
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We introduce a one-dimensional model basing on a special kind of asymmetrical one-side Ising-like 
dynamics. The model is equipped with the growth feature by adding a new node to the chain in each 
time-step. A spin of every succeeding node is influenced by the previous one but not vice- versa (the 
system in not hamiltonian). Owing to model's simplicity (it is an equivalent to a Markov chain) we 
are able exactly solve it in the presence of external field h. We conceive it as suitable for modeling 
for emotional online discussions arranged in a chronological order, where spin in every node conveys 
emotional valence of a subsequent post. We study the dependence of the chain's average spin (s) on 
external influence h (community tendency toward a selected valence) and noise level T (uncertainty 
of the emotional behavior). This leads us to an observation of three distinct phases - the first, where 
discussion evolution is determined by its launching message only, the second, in which the mostly 
observed emotion is coping the external influence, a finally, the third one, where the outcome is 
subject to fluctuations. 

PACS numbers: 89.20.Hh, 89.75.Hc, 89.75.Da 



Due to its simplicity and fully analytical treatment, 
one-dimensional models serve physicists as a very useful 
and comprehensible tool. Of the exceptional importance 
backed by feasibility of calculations is the Ising model 
Its basic version that considered only nearest neigh- 
bors coupling has been recursively modified to include, 
for example, long-range interactions Jij ~ \i — j\~ a (i.e., 
no phase transition for a > 2 [2jversus the existence of 
phase transition for 1 < a < 2 0-l5(). One-dimensional 
Ising model has also been subject to analysis with respect 
to a random field [f| 0] , where it has been proved that the 
ratio of up spins to down spins exhibits devil's staircase 
Q. Other research touches the issue of domain statis- 
tics @, cluster size distribution [lj| or inhomogeneity in 
the external field [ill, EH- I n this paper we introduce 
a one-dimensional model basing on a asymmetrical one- 
side dynamics. However, the asymmetry is unlike the one 
proposed by Huang [l3l Il4{. where the spin variable can 
take on two eigenvalues +1 and — 1/A with A > 1 nor it 
is connected to degeneration of higher-energy spin state 
[loT ]. Instead, we modify explicitly Ising hamiltonian by 
taking into account only left-hand neighbor. Moreover 
we also equip our model with growing (or evolving) com- 
ponent, i.e., we grow the chain by adding one node in 
each time-step. 

Although extremely useful, one-dimensional models of- 
ten suffer from exaggerated simplicity, i.g., one finds no 
evidence to support the idea that the agents related to so- 
cial dynamics are to be distributed on a chain. In fact, in 
most cases, the relations between them tend to be highly 
embedded in a complex network space fl6j — reflecting 
the actual topology of human contacts (private or profes- 
sional). In this Letter we give clear reasons for choo sing 
this very topology. Motivated by recent research [I?], Hq 
we model the process of postings of emotional messages 
[l9l l20j in a Internet forum. It is a common situation 
that the participants of such media use a chronological 
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FIG. 1: (Color online). A scheme of the simulation process, 
(a-d) Consecutive steps of an exemplary simulation: (a) start- 
ing from an empty node, (b) adding random spin to the first 
node, (c) adding the next node, (d) inserting spin according 
to dynamics rule, (e) An effect of the simulation. Orange 
(gray, s = +1) and blue (dark gray, s — + — ) discs symbolize 
spins. 



structure of the incoming posts that can be easily re- 
garded as a one-dimensional chain (i.e., the consecutive 
posts are represented by the nodes in the chain). The 
results of our previous analysis [13, EH indicate that one 
of the most dominant phenomena seen in the online dis- 
cussions is the emotional homophily. It follows that the 
discussion's participants take into account the history of 
the thread and basing on that knowledge issue their own 
posts, trying to imitate the overall emotional tendency. 
The simplest possible way of implementing this behavior 
is to take into account only the last comment (i.e., the 
newest one). 

Model description. The model bases on the idea of 
a growing chain (see Fig. [IJ. The process is orga- 
nized as follows: in the first step, a random variable 
So = ±1, called spin (an equivalent of emotional va- 
lence [2lJ in online discussion), drawn with probability 
Pr(so = ±1) = 1/2 is inserted into the first node of 
the chain (Fig. QJi-b). Then, another node of the chain 
is added to the right of the last one (Fig. [TJ;) and it 
is initially equipped spin once again drawn with equal 
probabilities Pr(si = ±1) = 1/2. Subsequently the node 
becomes a subject to the updating procedure that is bas- 
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FIG. 2: (Color online) Average spin (s) as function of tem- 
perature T for different values of the external magnetic field 
h = 0.01 (squares), h = 0.1 (circles), h = 0.2 (upward tri- 
angles), h — 0.5 (downward triangles), h — 1.0 (diamonds), 
h = 2.0 (hexagons). Solid lines come from Eqs (((J and Q. 
All data points are for N — 10 3 and averaged over M = 10° 
simulations. 



FIG. 3: (Color online) Average spin (s) as function of tem- 
perature T for different values of the chain size N = 10 2 
(squares), 10 3 (circles), 10 4 (upward triangles), 10 5 (down- 
ward triangles), 10 6 (diamonds). All data points for ft = 0.1 
and are averaged over M = 10 5 simulations. Solid lines come 
from Eq. (|6]) and the dotted line is tanh(2/i/T). 



ing on the Ising-like model approach. For each new node 
i we define Si = — Je^ie,; — /ie,;, where J is an equivalent 
of exchange energy and h is the external field. The func- 
tion Si promotes spins of the same sign in the consecutive 
nodes of the chain and can be treated as an emotional 
discomfort function. As the spin is drawn, we test how 
flipping its sign to the opposite one (i.e., from s; = +1 
to Si = — 1 or likewise) affects the change of function S 
as AS = Si — Si = — (Jsi-i + /i)(s^ — s^, where term 
Si corresponds to Sj calculated when Si —> = — Sj. 
Then we follow the Metropolis algorithm [22| i.e., if the 
AS < we accept the change, otherwise we test if the ex- 
pression exp(— f3AS ) is smaller or larger than a random 
value £ £ [0; 1] (here (3 = -jzjpp, where fcs is Boltzmann 
constant and T is the temperature). If the latter oc- 
curs we accept the change otherwise the spin is kept as 
originally chosen. The procedure of adding new nodes 
and setting their spin variables according to the above 
described rules is repeated until size TV of the chain is 
reached (Fig. [TJ:). As an outcome, the average spin in 
the chain (an equivalent of the average emotion in the 
discussion), defined as (s) = YaLi s » i s calculated. It 
is then averaged over M realizations (typically, in this 
study M = 10 5 ). 

Numerical simulations. — All numerical simulations 
have been performed for J = fcg = 1- Figure [2] shows 
the average spin (s) as a function of the temperature T 
for selected values of external field h. In the case of h < 1 
the plot reveals (s) equal to zero for small T«l, then 
a clear maximum for some specific value T c appears. Fi- 
nally, a recoil toward zero for T > T c takes place. In 
the case of h > 1 such a phenomenon is not observed: 
instead, (s) = 1 for small T and then there is a decrease 
toward zero. 

The explanation of such behavior is following: if the 
field h is smaller than the value J = 1 of influence from 
the previous post and the noise level (temperature) T is 
low the (s) follows the spin value of the first node of the 



chain. Let us note that in our system the initial condi- 
tion for the node sq is also a system boundary condition. 
Taking into account the fact that there is an equal chance 
of starting from either so = +1 or sq = — 1 this leads to 
a situations where half of the realizations are (s) = +1 
and half (s) = — 1. While temperature increases the fluc- 
tuations cancel the effect of influence of the launching 
spin (memory about system initial conditions) and more 
and more spins tend to align according to the sign of the 
external field leading to the observed peak. Since fluctu- 
ations destroy also the ordering impact of the field h in 
the whole system thus after crossing some characteristic 
temperature T c the average spin is lowered. Since the 
effect results from an interplay between initial conditions 
external field h and the fluctuations thus the temper- 
ature T c should be small for larger systems when even 
small fluctuations destroy the initial condition influence 
(it will discussed later). If the field h is stronger than 
the influence from the launching node (h > 1) then every 
next spin follows this field thus thermal fluctuations T 
can only lower the mean spin (s). Such a behavior is in 
fact observed at two upper plots of Fig. [5J 

Figure [3] shows that indeed for smaller systems (e.g., 
= 10 2 the characteristic temperature T c is of order 
0.1 — 1 and is shifted toward lower values as compared to 
systems of size N = 10 6 . An observed maximum (s)(T c ) 
is smaller for smaller systems since the influence of initial 
conditions plays for them a more important role. 

It is also interesting to track the dependence of the 
average spin value on the level of the external field (see 
Fig. 2]). In the case of low temperatures (T < T c ) average 
spin value changes abruptly from (s) = — 1 to (s) = for 
h = — 1 and then from (s) — to (s) = 1 for h = 1. For 
higher temperatures (T > T c ) this change is smoother 
and length range of h for which (s) = is smaller. 

Analytical description. — For modeling purposes we 
start with checking the way the probability of positive 
emotion following another positive emotion is calculated 
in the presence of the external impact h > 0. First, we set 
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while for h > J wc obtain 



FIG. 4: (Color online) Average spin (s) as function of the ex- 
ternal field h for different values of temperature T — 0.01 
(squares), T — 0.1 (circles), T = 0.2 (upward triangles), 
T = 0.5 (downward triangles), T = 1.0 (diamonds), T = 2.0 
(hexagons). Curves for (a) N = I0 2 , (b) TV = fO 5 . All data 
points are averaged over M = 10° simulations. Solid lines 
come from Eq. ((6]). 




FIG. 5: (Color online) Characteristic temperature T c versus 
the size of the chain N. Symbols (upward triangles for h = 
0.01, downward triangles for h = 0.1, diamonds for h — 0.5 
and circles for h = 0.9) are numerical solution of Eq. (j6j) 
while the solid line comes from Eq. Q. 



s = 1. Then, with equal probabilities Pr(si ± 1) = 1/2, 
the spin in the next node is chosen to be positive or neg- 
ative. Next, we calculate the change of function S : (i) if 
si = +1 and si = -1 then AS = 2(h + J) > 0, so the 

change is accepted with probability equal to e~P(. h+J J 
and not accepted with probability equal to 1 — e~^ h+J ', 
where p = 2/(k B T), (ii) if s 1 = -1 and s[ = +1 then 
AS = —2(h + J) < 0, so the change is always ac- 
cepted. As a consequence the probability p ++ of pos- 
itive spin following another positive spin is equal to 



P++ 



1(1" 



-P(h+J) 



)+hxl = 1- 



-P(h+J) 



. Following 



the same line of thought it is possible to obtain the ex- 
act expressions for the probabilities p ++ ,j>_| ,p y,p 

depending on the range of h. Then, for h G [0; J) wc 
have 
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For simplicity of the further calculations we use the 
following notation: wc denote the probability to stay in 
+1 state as p++ = p, the probability to move from state 

+1 to state —1 is p^ = 1 — p. Similarly wc denote the 

probability to stay in — 1 as p = q and the probability 

to move from state —1 to +1 is p |_ = 1 — q. In effect we 

obtain the transition matrix P 



P 1 ~P 
1-q q 



(3) 



that defines probabilities evolution of both states 
Pr(s„ = ±1) as s n+ i = s„P, where s„ = 
[ Pr(s„ = +1) Pr(s„ = -1) ] . Thus the evolution of s„ 
is in fact an equivalent of a two-state Markov chain (23j 
governed by the transition matrix P that has a solution 
s„ = SqP", where so = [ 1/2 1/2 ] and 



pn 



g-l + (p-l)(g+p-l)" (p-l)[l-(g+p-l)" 



(9- 



q+p-2 
-l)[l-(g+p- 



q+p-2 



q+p-2 
-l + (g-l)(9+p-l)" 
q+p-2 



(4) 



This leads to calculation of the average spin in the n- 
th node as (s„) = s n [ 1 — 1 ] . Finally, performing the 
mean of (e„) over all nodes in the chain gives the average 
spin in the chain: 



(s) 



p-q 

-p-q 



1 



1 

N 



l-(p + q-l) N+1 



N(2-p-q) 



, (5) 



Similar calculations can be performed for ranges h G 
(— oo; — J] and h G ( — J;0]. The symmetry of the prob- 
lem results in swapping all the indices "+" to "-" and 
likewise in Eqs ([T][2]) . As an outcome we obtain a rotated 
matrix P that leads again to Eq. ([5]). In effect by apply- 
ing the of exact values of probabilities p and q given by 
Eqs ((T][2|) and obtain the relations that bind the average 
emotional value with external field and temperature for 
ranges \h\ < J: 



(s) s = tanh/3/i 
and \h\ > J: 



1+ N- 



l 1 - (l -e~ /3J cosh/i/i 



N+l 



Ne-' 3J cosh j3h 



(6) 



(s)i 



cosh P J — c 
sinh P J — c 



8h 



1+ N- 



-? h sinh PJ 



N+l 



Nil- e-P h sinh pj 



. (7) 

The above equations are subject to further approxima- 
tion. For example, (s) s w tanh/3ft, for large N (see dot- 
ted line in Fig. [3]), while for large both and T we have 
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(s) w f3h for the whole range of h. Equations <[6][7]) fully 
support (solid lines in Figs [2] and 0]) numerical observa- 
tions. Both — numerical and analytical — approaches 
indicate crucial role played by the system's size N (see 
Fig. [5j : for a constant value of external field h increas- 
ing N leads to a shift in T c toward T = as well as to 
increase of the value (s)(T c ). It is essential to consider 
its exact value; first of all, as it has been shown before 
in the previous Section this phenomenon is observed only 
for h < 1. Then the formal solution can be obtained from 
Eq. © but it is only possible to treat it numerically (see 
symbols in Fig. [SJ To get an analytical estimation we as- 
sume that (3h <C 1 which gives the opportunity to rewrite 
Eq. © as 



1 

TV 



1-1- e-? J 



N+l 



Ne~P J 



(8) 



The above expression is an equivalent of magnetic suscep- 
tibility (j^-Qj^J times h. If we further assume N ^> 1, 

and solve d g S ^ B = one can approximate T c as 



2 J 



k B [W(iVc) - 1] 



(9) 



where W(...) is Lambert W function. Comparison be- 
tween this approximation and numerical solution of Eq. 
© is shown in Fig. [51 providing evidence of good agree- 
ment for small values of h as expected. If follows that for 
real-world examples where the discussion length prob- 
ability follows a power-law function p(N) ~ TV~ 7 and 
N-max ~ 10 5 [13] the effect of non-zero T c could be ob- 
served. 

Discussion. — It is of use to compare and contrast the 
obtained results with the classical one-dimensional Ising 
model (e.g., [13 )■ The most noticeable difference regards 
undoubtedly the foundations of the model — in the clas- 
sical case the length of the chain is fixed and each node i 
is initially filled with a spin s; = ±1. The thermalization 
procedure involves energy couplings with both neighbors 
of node i and it is performed repeatedly. Here, in the 
presented model, the chain grows by adding new nodes 
and each node is subject to the thermalization procedure 
only once, moreover just its left-hand neighbor is taken 
into account. The magnetization per spin (an equivalent 
of (e)) in the case of the classical one-dimensional Ising 



model is given by 



m(h,T, N) 



sinh j3h 



l-f$= 



N 



^sinh 2 /3h + e~ A ? J 1 + 



N 



(10) 



( cosh/3/i ± yjc-^- 1 + sinh 2 /3h) arc the 



where A± 

eigenvalues of the transfer matrix [24|. In our model 
the size of the chain plays a pivotal role, governing the 
position and height of the maximum of e for h < J. 
Although Eq. ([9]) proves that in the limit of N — > oo 
the value T c —> 0, nonetheless the convergence is slow 
(approximately as 1/ In N) . For the classical Ising model 
the term containing system's size rapidly converges to 1: 
e.g., for J = ks = 1) h = 0.1 and T = 0.1 one needs as 
little as N = 10 to have the convergence. On the other 
hand, in our model a large chain length is needed. For 
small T the terms in square brackets in Eq. ([6]) can be 

approximated by \{N + l)e~^ J cosh f3h which in turn can 
be used to estimate the critical value of for which / 

is equal to one as N c « 4c^ J -' l \ Then, for the same 
parameters (J = k B = 1. h = 0.1. T = 0.1) N c = 2.62 x 
10 8 . 

Finally, opposite to our model, the magnetization in 
the classical Ising model for J > (ferromagnetic case) 
is a strictly monotonous decaying function of T starting 
from m = 1 for h > 0. The maximum is observed for 
J < (antifcrromagnctic case) but it is independent from 
TV (in the sense of factor's // rapid convergence to 1) and 
for ftclwc have T c = 2. 

Summary. — In this paper wc proposed and exactly 
solved a one-dimensional model basing on a growing 
chain with an asymmetrical Ising-like dynamics. The 
model can be considered as a certain modification of the 
original Ising concept accomplished by restriction im- 
posed on the hamiltonian function. Moreover, we also 
introduced growth factor into the model, expanding the 
chain with one node in each time-step. Those features en- 
abled simple analytical treatment (as Markov chain) that 
produces exact expressions for an equivalent of magneti- 
zation in the classical Ising model. The model, thought 
as a construct resembling on-line discussions, combines 
three different elements (i) chain-like structure of on-line 
discussions, (ii) tendency to take into account only the 
previous comment, (iii) global interactions inserted in a 
form of external field and fluctuations. 
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